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Abstract 

Elastic waves with transverse polarization inside incidence plane can create longitudinal surface 
wave (LSW) after reflection from a free surface. At a critical incidence angle this LSW accumulates 
energy density, which can be orders of magnitude higher than energy density of the incident 
transverse wave. A specially arranged vessel for storage of ultracold neutrons (UCN) can be 
used to verify this effect. 
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I. INTRODUCTION 



It is knowni- that reflection from an interface of an elastic wave of a given mode in 
isotropic and anisotropic media is accompanied with mode conversion. Inl it was noticed that 
at grazing angles smaller than some critical angle the mode conversion can result in creation 
of a longitudinal surface wave (LSW), which is alike to those appearing at earthquakes. The 
most prominent property of LSW is that critical angle it can accumulate an energy density 
orders of magnitude larger that that of the incident plain wave. 

In the next section the theory of elastic waves in isotropic media is presented according 
to^, and in the third section an experiment with ultracold neutrons (UCN) (see, for example^) 
is proposed to observe effect of LSW on UCN storage. 



II. ELASTIC WAVES IN ISOTROPIC MEDIA 

The main element of the theory of elastic waves 9 is a displacement vector u{r, t) of 
a material point at a position r at a time t. Its Cartesian coordinates Ui(r,t) obey the 
Newtonian equation of motion 

d 2 d 

where p is the mass density, Xj are components of the position vector r, is a stress tensor 

a ij = c ijklUkh (2) 

which is proportional to the deformation or strain tensor Uij 

1 / duj du k \ 

and coefficients of proportionality Cijki in fl2]) form a tensor which in isotropic media looks 

Cijki = XSijS k i + n{5u5 k j + S ik 5ij), (4) 

where 5ij is the Kronecker delta function, and A and // are Lame elastic constants. 

With the stress tensor (J2J), (Jl|) the Newtonian equation of motion ([ID for the displacement 
vector becomes 

P Ui = VjCTij = p[A Ui + Vi(V • u)] + AVi(V • u). (5) 



Usually the displacement vector u{r, t) is represented as a sum of two parts u(r, t) = 
V(/? + V x ip, where (p is a scalar and tp is a vector potentials. This tradition was broken 
in -, and it helped to find new unexpected effects. The solution of §5§ is represented in the 
form of a complex plane wave 

u(r, t) = u A exp(ikr — iut), (6) 

where vector A ia a unit polarization vector, and Uq is the wave amplitude with dimension 
of length. After substitution of u from ([6]) into (J5j) one obtains an equation for A: 

pu 2 A = pk 2 A + {\ + fj,)k(k- A). (7) 

Because of the linearity of the equation the amplitude uq of ([6]) does not matter and for 
some time will be omitted. 

For a given propagation direction k = k/k one can introduce two orthonormal vectors 
and e^ 2 \ which are perpendicular to k, and the orthonormal basis e^, e^ 2 \ = n, 
in which the polarization vector A looks 

A = a^e^+a^e^ + a^e^. (8) 

After multiplication of (JTj) by each of the one obtains three equations for the coordinates: 

(pu 2 - fik 2 )a {1 ' 2) = 0, (pu 2 - [A + 2fi]k 2 )a^ = 0. (9) 

They are independent and give three solutions: two for shear modes A 1,2 = e^ 1 ' 2 ) with speed 
Q = ' p and wave number k^ 1,2 ^ = u/ct, and one longitudinal mode A 3 = k, with speed 
c\ = J(A + 2fi)/p and wave number k 3 = u/q. 

Consider now reflection from a free surface at z = of the wave (J6j) propagating in the 
medium at z < 0. In this case it is convenient to choose basis vectors e^ 1,2 ) _L n in such a 
way, that is perpendicular to the incidence plain, and lies in it. The most interesting 
is reflection of the shear mode A^ = e^ 2 \ The displacement vector, when the incident wave 
is of this mode, looks 

u(r,t) = exp(ife||7*|| — iut) 

where r^ 2 ' are reflection amplitudes of the mode 2 with transformation into mode j = 
2,3, and kj_* = J k^ 2 — fc|. It is important to notice that k^ < h±\ therefore reflected 
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A (2) eifc f z + r ( 22 ) A g) e _ jfc f, 



+ r^Afe- ik 



(10) 



longitudinal wave propagates in nonspecular direction. Its grazing angle cp^ is less than ip^ 
of the incident wave. Therefore at some critical angle tp^ = arccos(ct/cz) the longitudinal 
wave starts to propagate along the surface, and at ip^ < (p^ the component k^ becomes 
imaginary, i.e. the longitudinal reflected wave becomes LSW. This is the most interesting 
phenomenon which, is a model of the earthquake waves, because the most devastating effect 
of earthquakes results from vibrations along the surface. 

Let's find the energy of the LSW 10 . To do that we need to know reflection amplitudes. 
They are found from the boundary condition, which requires the stress vector T with com- 
ponents Tj = (TjiTii, where rij are components of a unit vector n along the normal, to be 
zero at z = 0. For stress tensor (T5]), (j4]) the stress vector for the displacement u(r,t) is 

T(u(r, tfj = An(V ■ u) + fj,[V(u ■ n) + (n ■ V)w]. (11) 

Substitution of ( TTOj) in it gives boundary condition in the form 

s (2) +r (22) s (2) +r (32 )£ j(3) =0) (12) 

where for every plane wave Aexp(ikr) the vector B is defined as 

B = Xn(k- A) + fi[k(A • n) + A(n • k)], (13) 

is related to the incident wave with k^ = ku + nkj_ , and B^ (j = 2,3) are related 
to the reflected waves with k^ = feii — nk^ . 

Multiplication of ( !T2l) by two unit vectors: n and r = fe||/|fe||| gives a system of two 
equations for determination of r^ 2 \ Their solution is 

r (32) _ ^ 4fcffc||(fc( 2 ) 2 -2fcg) r(22) _ Akfkfkl - {k^ - 2k\f 

" Akfkfhl + (A;( 2 ) 2 - 2A; 2 ) 2 ' ' ~ Akfkfk 2 + (A;( 2 ) 2 - 2fc 2 ) 2 ' 

and it is seen that imaginary k^ makes \r <y22 * > \ = 1, which means the total reflection of the 
mode 2. 

Because of energy conservation the energy flux density of the incident wave along the 
normal to the interface must be equal to the sum of energy flux densities of the reflected 
waves. The energy flux is a real quantity, and with our complex solutions it is defined 10 as 
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*dui du\ 



Re 



'"It 



(15) 
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where Re(F) means real part of F, and * means complex conjugation. After substitution of 
all the quantities we obtain 



J = Re(puciu 2 [EA*(k* ■ A) + k*}), 



(16) 



where E = 1 + A//i, and the amplitude uq of the oscillations is taken into account in fl6]). 
Expression ( fl6|) is valid for real and complex wave vectors and polarizations. The absolute 
value of this flux for the incident transverse wave is 



J = ptuc t u k. 

The flux of the LSW along the surface with account of its factor r^ 32 ) is 

= pum 2 |r (32) | 2 c 2 fe|| exp(2K lZ ). 



The ratio Q = Jg (z = 0)/J of two fluxes after substitution of r^ 32 ) from ( TH1) is 



or 
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fen 16fe[ 2) 2 £f(fe( 2 ) 2 - 2fe 2 ) 2 
k^WKfk^ 2 k\ + (fe( 2 ) 2 - 2fe 2 ) 4 ' 

cos ip sin 2 (4(/j) 
4(cos 2 if — cos 2 ip c ) sm 2 (2ip) cos 2 (<p) + cos 4 (2(£>) ' 



fen r 2 
n _ K \\ C l I (23)|2 

W fe( 2 ) C 2 ' ' 







(17) 



(18) 



(19) 



(20) 




FIG. 1: Dependence of Q(x) on x = cos if calculated for E ~ 2.3 like in stainless steel. The left 
point on the abscise axis corresponds to critical cos(/? c = q/q ~ 0.57, where Q = Q c = 16. In the 
case of E = 1.2 the similar calculations give cos</? c = q/q = 0.674 and Q c = 323.616. 
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where if = {p( 2 \ Dependence of this function on cos<£> is shown in Fig. [TJ We see that the 
highest energy density is accumulated in longitudinal surface wave, when (p is slightly less 
than ip c . The higher is the ratio q/q, the smaller is the critical grazing angle and the higher 
is accumulation of the LSW energy density near the critical angle. In the case of stainless 
steel the ratio q/q is near 0.56, and Q(<p c ) ~ 16. 

a. Critical polarization. Let's nevertheless find the exact polarization of the total wave 
field on the interface, when the grazing angle of the incident wave is ip c . According to ( II Op 
the total displacement vector at z = is equal to 

A t = A< 2 > + r^\ Vc )A^ + r^\ Vc )A%\ (21) 

where according to ( Fl4l) 

r^ 2 \^ c ) = 2sin(y? c )tan(2^ c ), r^ 22 \ Vc ) = -1. (22) 
Substitution of ( 1221) and 

A (2) = Tcos(y? c ) — nsin(<£> c ), A$ = — [rcos(^ c ) + nsm((p c )], Aft = r (23) 
into ( 12T1) gives 

A t = 2 cos((p c ) [t&n((p c ) tan(2v? c ) - l]r, (24) 

where sin((^ c ) = q/q. 

From this expression it is seen, that at critical angle vibrations are along the surface 
(such vibrations are the most devastating at earthquakes), and their amplitude becomes 
especially large at (p c w 7r/4, i.e. in the case when q/q ~ 0.7. In the case of stainless steel 
with Ct/ci = 0.57, though Q(ip c ) = 16 as is shown in Fig. CD the total amplitude is only 
A t (<f c ) — l-4r. In the case of q/q = 0.65 we get A t ({p c ) = 6.8r and Q((p c ) = 71 

At tp < Cf/ci polarization of LSW is directed along the wave vector fc*- 3 -* = (k\\,iK±), but 
this vector has an imaginary normal component K^. It means that vibrations in LSW have 
also normal component, which is shifted by phase it/ 2 with respect to vibrations along the 
surface. However near the critical angle the component K± is small and vibrations along 
the normal can be neglected. 
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III. AN EXPERIMENT WITH UCN TO OBSERVE LSW 



Scheme of the proposed experiment with UCN is shown in Fig. [2j The floor of the 
storage box for UCN is a thick basement. Several vertical thin plates are attached to it. 
And transducers of shear waves with polarization in the vertical plane generate ultrasound 
waves going to the storage floor at different angles. The UCN can be stored in the box or 



storage box 



UCN 




transducers 

FIG. 2: Scheme of an experiment with UCN sensitive to LSW. 

continuously flow through the box, and dependence of the storage time or transmission of 
the box on angle of incidence of the ultrasound wave is measured. If the grazing angle of the 
incident ultrasound wave is close and slightly below the critical one, the LSW created on 
the surface will make oscillating the vertical plates, and the UCN in the box will be heated 
at every collision with the plates and continuously increase their energy up to the limit after 
which they will go away through the box walls. It decreases the storage time or transmission 
of the box. 

For estimation of the effect one can consider some material like stainless steel with density 
p ~ 8 g/cm 3 , the speed of transverse waves c t ~ 3 km/s and c t /q = 0.57, like that used in 
Fig. [TJ Then ultrasound with cu — 1 MHz and amplitude u = 0.1 /i, has energy density 
pco 2 cu 2 of the the order of 10 W/cm 2 . The speed vq of the horizontal vibrations on the box 
floor surface at critical incidence angle will be near v\ = 0.4 m/s. It will create vibrations 
of the vertical plates with the same speed. Every collision of UCN with walls will increase 
in average their v 2 by vf = 0.16 m 2 /s 2 . Therefore, after 100 collisions with the walls their 
energy will be high enough to penetrate the walls and to escape from the storage box. 

The dependence of the storage time r or transmission of the box on cosine of the grazing 
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angle ip is expected to look as shown in fig. |3j Observation of the dip will prove the 
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FIG. 3: Qualitative dependence of the storage time t on cos(/j, where <p is the grazing angle of the 
incident sound wave. The position of the dip depends on ratio of q/q of the material used for 
basement in Fig. (2j The depth and the width of the dip depend on the power and frequency of 
the ultrasound. 

appearance of the LSW and therefore will prove the validity of the earthquake model. 

IV. CONCLUSION 

The theory of LSW is presented. Though surface waves were studied extensively^ the 
longitudinal one was discovered 7 only because of nonstandard approach to theory of elas- 
ticity. The LSW appears at reflection from a free surface of a shear wave with polarization 
in the incidence plane, when grazing angle of the incident wave is below of the critical one 
tp c = arccos(c t /c/). Below and near ip c the energy density of LSW can be orders of magni- 
tude higher than energy density of the incident wave. From the properties of LSW one can 
conclude that they can have the same devastating effect as earthquakes. Therefore LSW 
can be considered as a model of earthquake waves. We described here an experiment, which 
can check the properties of LSW. It is related to storage or transmission of UCN. In the 
experiment one can see a sharp drop of storage time or transmission of UCN through the 
storage box when LSW appears on the surface of the storage box floor. 
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V. HISTORY OF SUBMISSIONS AND REJECTIONS 



The paper was submitted to Phys.Rev.Lett. on 07.11. On 10.11 I received a short letter 
from Senior Assistant editor with words that my paper is not suitable for PRL. I submitted 
immediately to Europhys.Lett and and this time the fate of the paper was extraordinarily 
lucky. On 07.12 I received a positive referee report, and the editor (it was Prof. Rudolf 
Treumann) returned me my manuscript with completely corrected English. I cannot with- 
stand the desire to express my admiration and deep gratitude to Prof. Treuman. 

In the new version, submitted to arXiv I added a paragraph at the end of section II. 
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